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ABSTRACT 

We have used the Robert C. Byrd Green Bank Telescope to time nine previ- 
ously known pulsars without published timing solutions in the globular clusters 
M62, NGC 6544, and NGC 6624. We have full timing solutions that measure 
the spin, astrometric, and (where applicable) binary parameters for six of these 
pulsars. The remaining three pulsars (reported here for the first time) were not 
detected enough to establish solutions. We also report our timing solutions for 
five pulsars with previously published solutions, and find good agreement with 
past authors, except for PSR J1701— 3006B in M62. Gas in this system is proba- 
bly responsible for the discrepancy in orbital parameters, and we have been able 
to measure a change in the orbital period over the course of our observations. 
Among the pulsars with new solutions we find several binary pulsars with very 
low mass companions (members of the so-called "black widow" class) and we are 
able to place constraints on the mass-to-light ratio in two clusters. We confirm 
that one of the pulsars in NGC 6624 is indeed a member of the rare class of non- 
recycled pulsars found in globular clusters. We also have measured the orbital 
precession and Shapiro delay for a relativistic binary in NGC 6544. If we assume 
that the orbital precession can be described entirely by general relativity, which 
is likely, we are able to measure the total system mass (2.57190(73) Mq) and 
companion mass (1.2064(20) Mq), from which we derive the orbital inclination 
(sini = 0.9956(14)) and the pulsar mass (1.3655(21) M0), the most precise such 
measurement ever obtained for a millisecond pulsar. The companion is the most 
massive known around a fully recycled pulsar. 
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Subject headings: globular clusters: individual (M62, NGC 6544, NGC 6624)— 
pulsars: individual (J1701-3006D, J1701-3006E, J1701-3006F, J1807-2459A, 
J187-2500B, J1823-3021C, J1823-3021D, J1823-3021E, J1823-3021F) 



Introduction 



Millisecond pulsars (MSPs) form by "recycling" a dormant neutro n star through th e 
accretion of mass and angular momentum from a binary companion (lAlpar et al.l Il982l ). 
This leads to a very rapidly rotating and stable pulsar with a relatively low magnetic field 
(~ 10^ G) and long lifetime (> 10^ yr). MSPs form naturally in the dense environments 
of globular clusters (GCs), thanks to fre quent exchange intera ctions that may lead to the 
formation of mass transferring binaries ( Camilo fc Rasidl2005l ). Sensitive searches of GCs 
have uncovered 144 pulsars in 28 clusterqj, and the vast majority of these are recycled MSPs. 
Indeed, nearly half of all MSPs have been discovered in 

The same interactions that form MSPs so efficiently in clusters also lead to many exotic 
systen as that are rarely s een in the disk of the Gal axy. These include the fastest spinnin g 
MSP ( iHessels et al.ll2006), highly eccentric binaries ( iRansom et al.ll2005l : iFreire et al.l 120071) 



massiv e neutron stars (IFreire et al.l 120081 ) , pulsar-main se quence binaries ( iD'Amico et al. 



2001bl ). and many "black widow" systems (iKing et al.ll2005l ). These discoveries demonstrate 
the huge scientific payoff that can come from the discovery of unique pulsars thanks to two 
factors — the extreme nature of neutron stars, which opens windows on otherwise inaccessible 
realms of physics, and the extraordinary clock-like nature of MSPs. The bedrock of pulsar 
astronomy is timing, the process of creating a model that unambiguously accounts for every 
rotation of the pulsar, thus probing the pulsar and its environment. For MSPs, the arrival 
time of a pulse can typically be measured to within a few microseconds or better (e.g. 
Verbiest et al.ll2008l ). which enables very precise timing models. Timing GC pulsars leads to 
unique challenges, though. Because GCs are usually at a distance of several kiloparsecs, the 
fiux density of their constituent MSPs is usually very low, necessitating very long integration 
times. Another consequence of cluster distances is very high dispersion measures, which 
necessitate moving to higher observing frequencies (e.g. 2 GHz) where pulsars tend to be 
weaker. Luckily, since many clusters contain several MSPs that can each be observed during 
a single observation, the required observing time is well spent. 



See http://www.naic.edu/~pfreire/GCpsr.htiiil for an up-to-date list. 

^There are 160 field pulsars with period P < 20 ms in the ATNF catalog, and 131 GC pulsars that meet 
the same criteria (see http : //www. atnf . csiro . au/people/pulsar/psrcat/| ). 



- 3 - 



In addition to recycled MSPs, a s mall population of slow, non- recycled pulsars have also 
been observed in a handful of clusters (ILyne et al.lll996l : iBoyles et al.ll201l[ ). Their presence 
is something of a mystery, because they resemble in every way the "normal" pulsars that 
are so numerous in the Galactic disk, but which have lifetimes ~ 10''-10^ yr and are thought 
to form through core collapse su pernovae of massive stars. However, GCs are old stellar 
systems with typical ages 10^° yr ( ICarretta et al.l 120001 ). so all stars massive enough to form 
pulsars should have died some 10^° yr ago. As such, the core collapse of massive stars 
cannot be the avenue through which non-recycled pulsars in GCs have formed. The most 
popular alternati ve formation s cenar io involves the collapse of a massive white dwarf via 
electron capture ( lNomotdll984j . 119871 ). though the details of these so-called electron capture 
supernovae are not well understood. 

In this paper, we present new timing solutions for nine pulsars spread across three 
GCs— M62, NGC 6544, and NGC 6624. Each of these clusters contains multiple MSPs and 
NGC 6624 contains two non-recycled pulsars. All of these pulsars were discovered elsewhere, 
but full, phase-connected timing solutions have only been published for five of them (see 
Table [T] for a summary and list of references). We also present our solutions for these 
five. We began a timing campaign using the Robert C. Byrd Green Bank Telescope (GBT) 
for the nine pulsars without full solutions, starting in 2009 February. We have managed 
to obtain full solutions that include measurement of the first period derivative for six of 
these pulsars. A combination of low flux densities and an irreversible data processing error 
conspired to prevent us from obtaining full solutions for the two remaining MSPs, but we are 
in a position to comment on them in further detail. Our most exciting result is the discovery 
of a massive companion orbiting a fully recycled MSP. In ^ we describe our observations 
and in ^ outline our method for timing the pulsars and further data analysis. A discussion 
of individual systems can be found in ^ and we provide a summary in ^ 



2. Observing Scheme 



Observations were carried out with the GBT, observing at a frequency of 2.0 GHz using 
800 MHz of bandwidth, although persistent radio frequency interference (RFI) reduced the 
usable bandwidth to ~ 600 MHz. This frequency was chosen to overcome the deleterious 
effects of dispersive smearing caused by free electrons in the ISM, and has been used success- 



fuUy by our group before. Data were recorded in the PSR FITa'^l (iHotan et al.l 20041 ) format 
using the Green Bank Ultimate Pulsar Processor (GUPPI) (jPuPlain et al.ll2008l ) with a 64 /xs 



http : //www. atnf . csiro . au/reseaxch/pulseir/psrf its/index. html 
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sampling time and 2048 frequency channels across the entire bandwidth. We typically ob- 
served NGC 6544 for 30 minutes, and M62 and NGC 6624 for 45-60 minutes each, although 
the exact integration times varied between observations. In general, data were relatively free 
of RFI, but when necessary, we used the RFI excision tools in the PRESTC0 software suite to 
mask out contaminated portions of the data. 

We also used ar chival GBT data ta ken in 2004. These data were collected using the 



GBT Pulsar Spigot flKaplan et all 120051 ) at either 820 MHz (with 50 MHz of bandwidth) 
or at 2 GHz (with the same bandwidth as above). Most of these observations used 2048 
frequency channels and 40.96 /is sampling, although some used 1024 channels and 81.92 /is 
sampling. All of these data were analyzed using PRESTO. 

The pulsars in NGC 6544 were observed again in 2011 as part of a campaign to measure 
Shapiro delay in PSR J1807— 2500B (see 14.61) . These data were collected using GUPPI in 
a coherent de-dispersion search mode (i.e. a filterbank where each channel was coherently 
de-dispersed) at 1.4 GHz using 800 MHz of bandwidth. For these observations we used 512 
frequency channels and 10.24 /is sampling. 



3. Timing and Data Analysis 

Data were processed using a combination of PRESTO and PSRCHIVE0. All of the folded 
pulse profiles were phase-aligned and summed to create high signal-to-noise (S/N) profiles. 
The profiles were fit with one or more Gaussians, from which standard pulse templates were 
made and used to obtain pulse times of arrival (TOAs) via cross-correlation in the Fourier 
domain. We typically obtained three to six TOAs per pulsar per observation depending on 
the S/N of that observation, though some pulsars (such as NGC 6544A) were sufficiently 
bright that many precise TOAs could be measured. Tim ing solutions were co nstructed 



by performing a weighted fit to the data using TEMPOaill f lEdwards et al.l 120061 ) with the 
DE405 Solar System ephemeris and TT(BIPM2011) time standard. All values are reported 
in Barycentric Dynamical Time. It is not uncommon for timing models to have a reduced 
> 1 (Xred) even after all parameters have been well measured. When no systematic trends 
are present in the data, it is assumed that this value of Xrcd due to an underestimate of the 
error on individual TOAs. As is common practice, we deal with this situation by multiplying 



' http : //www . cv . nrao . edu/~srEinsoin/presto/ 



http : // psr chive . sour cef orge . net /| 



If 



http: / / sourcef orge . net/pro jects/tempo2/| 
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TOA errors by a small constant error factor so that Xred ~ ^■ 



3.1. Corrections for Cluster Acceleration 



MSPs have a small intrinsic rate of spin-down (Pint) that is usually heavil y contaminated 
by accelerations within the potential of the cluster and Galaxy ( lPhinneylll993l ). This makes it 
impossible to measure their spin-down related properties (surface magnetic field, spin-down 
luminosity, and characteristic age) directly. Instead, we calculate the limit 



int ^ -fobs ^ 0,c,max ^ '^G ^ A'" 



P 



P 



where P is the pulsar period, a^max and are the accelerations due to the cluster and 
Galaxy, respectively, fj, is the proper motion, D is the distance to the cluster, a n d c is the 
speed of ligh t (the last term accounts for the Shklovskii effect ( IShklovskiil 1 19701 ) ) . iPhinney 
( I1992I . Il993[ l showed that to within 10% accuracy. 



3al 



2c(i?2 + i?2^Jl/2 



(2) 



for Rpsr < 2-Rc, where is the central velocity dispersion, Rc is the core radius, and R 



psr 



is the projected distance of the pulsar from the center o f the cluster . Beca use all three of 



the clusters studied here are core collapsed, analytical (IFreire et al. 



2005) and numerical 



models are not applicable, which is why we use the approximations of 
developed fo r use with pu lsars in M15, which is also core collapsed. We use the va^ 
Rr fou nd in HarrisI ( 1996 . 20 10 edition ). The following references are used for cr,,: 



Phinneyl . These were 
ues of 



Harris 



(119961, 2010 edition) for M62; IWebbinkI (ll985|) for NGC 6544; and IValenti et all (120111 ) for 
NGC 6624. Table [2] lists the relevant properties of each cluster. The Galactic contribution is 
calcul ated under the approximation of a spherically symmetric Galaxy with a flat rotation 
curve (jPhinneylll993l ) and is 



Gal 



-7 X 10-^^ ( - 



COS b cos i + 



6 — cos b cos i 



1 + 5^ — 2(5 cos b cos i 



(3) 



where b and £ are the Galactic latitude and longitude of the cluster, 6 = Rq/D and Rq is 
the Sun's Galactocentric distance. However, the cluster term is usually dominant. 
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3.2. Flux Calibration and Rotation Measure 

Rough mean flux density (5*^) estimates were made by assuming tliat tlie off-pulse RMS 
noise level was described by the radiometer equation, 



a = , (4) 

GA/ripoiAz/tobs 

where Ttot is the total system temperature, G is the telescope gain, n^o\ = 2 is the num- 
ber of summed polarizations, and Az/ is the bandwidth. For the GBT 2 GHz receiver, 
G = 1.9 K Jy~^ and T^ys ~ 23 K + Tgky, where Tsky is the contributi on from the Galactic 



synchrotron emission. This was calculated by scaling the values from iHaslam et al.l (119821 ) 
with a spectral index of —2.6. The typical uncertainty in these estimates of Sy is 10%-20%. 
We were able to record full polarization data for one observation per cluster. These were 
calibrated using the 25 Hz noise diode on the GBT, and we searched for a significant rota- 
tion measure (RM) by looking for a peak in the polarized flux. We searched from —1000 
to 1000 rad m~^, but in most cases, the S/N was not sufficient to constrain RM. The no- 
table exceptions are the two pulsars in NGC 6544, for which we measure an average RM of 
~ 158 radm~^. Figure [1] shows the fully calibrated, RM corrected profiles of NGC 6544A 
and B. 



4. Discussion of Individual Systems 

All of our timing solutions can be found in Tables |3H71 along with some derived prop- 
erties of the pulsars. Average pulse profiles and Doppler modulated pulse periods of the 
binary pulsars in these clusters are shown in Figs. [2] & |3l Post-fit timing residuals are shown 
in Fig. m We discuss each individual system below. 



4.1. Pulsars with Previously Published Timing Solutions 

PSRs J1701-3006A, B, and C and J1823-3021A and B all have timing solutions pub- 



lished by other authors (jPossenti et al.ll2003l : iBiggs et al.l Il994j ). We constructed our own 



timing solutions based solely on our data to confirm that there were no irregularities in our 
data or methods. We find that nearly all our measured parameters agree with those pub- 
lished elsewhere to within errorj^. The one major exception to this is PSR J1701— 3006B 



^^Our measurement of right ascension for M62A differs from the results given in Possenti et aL ( 20031) by 
0.003" , which is approximately five times the formal 1-a errors 
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(hereafter M62B), where we see a highly sig nificant change i n orb ital period and dispersion 
measure (DM) compared with the results of iPossenti et al.l (120031 ) . We have also measured 
the rate of change of the orbital period Pb = —5.51(62) x 10~^^. We performed an F-test to 
determine if the addition of Pb was indeed required by the data. Without it, = 142.82 
with 71 degrees of freedom, while after fitting for Pb the improved to 73.94 with 70 de- 
grees of freedom. The probability that this improvement is due to chance is < 1.3 x 10~^^, 
so the measured Pb does indeed seem to be required by the dat a. M 62B is known t o eclipse 
and has an optical and X-ray counterpart (ICocozza et al.ll2008l ). As ICocozza et al.l discuss, 
the pulsar and companion star are almost certainly interacting. The contribution to Pb 
expected from general relativity (GR) is two orders of magnitude smaller than observed, 
assuming Mp = 1.4 Mq and Mc = Mc^min, and there is no realistic combination of pulsar and 
companion masses and inclination angles that would lead to such a large relativistic Pb. As 
such, classic tidal effects of the extended companion are probably the cause of the change in 
orbital period. Our R MS timing residuals are over a factor of two smaller than obtained by 



Possenti et al.l (120031 ) . which probably explains why we were able to detect Pb even though 



we had a shorter timing baseline. 



4.2. PSR J1701-3006D 



PSR J1701-3006D (hereafter M62D) is a 3.42 ms binary MSP. M62D (alon g with E and 
F) were discovered, and initial orbital solutions were given, by IChandlerl (120031 ). The orbital 



period is 1.1 days, with a small eccentricity I. e ~ 4.12 x 10~ . The minimum companion 
mass is ~ 0.12 Mq (assuming a 1.4 Mq MSP), and is likely a white dwarf. We see no evidence 
for eclipses in M62D, but none of our observations cover conjunction, when an eclipse would 
be most likely. Three observations do, however, start or end within four hours (~ 15% of the 
orbital period) of conjunction, and one observation ends only 15 minutes before conjunction. 
The lack of eclipses increase our confidence that the companion is a white dwarf, and not a 
main sequence (MS) star (for which eclipses should be common). We are unable to measure 
any precession in the longitude of periastron, so no further constraints can be placed on the 
mass or geometry of the system at this time. The accelerati on of M62D (PP~^) is somewhat 
higher than for the majority of GC pulsars ( lRansomll2008l ). though not exceedingly so. We 
note, though, that PSRs J1701— 3006E and F have similarly large accelerations. 



^•^W e have used the ELLl tuning model in TEMP02, which is appropriate when a sini/ce^ <^ 1 (jLange et al, 
200 ih . 
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4.3. PSR J1701-3006E 

PSR J1701-3006E (hereafter M62E) is a 3.23 ms binary MSP in a circular orbit. The 
orbital period is only 3.80 hr, and the projected semi-major axis is only 0.0302 R©. When 
combined with the low minimum companion mass (0.031 Mq, or o nly 31 Jupiter n iasses) 



and the presence of eclipses, M62E is clearly a black widow pulsar ( iKing et al.ll2005[ ). The 
presence of eclipses makes it likely that the system is being seen at a high inclination, so that 
the true companion mass is probably close to the minimum. We have good orbital coverage 
of the system, and observe both ingress and egress during eclipses, which are sharp and do 
not lead to a significant change in DM. The eclipses are centered around orbital conjunction 
and occur for ~ 12% of the orbit (or roughly 12 minutes). We see no evidence for eclipses 
at other orbital phases. The projected extent of the ecl ipsing mate rial is ~ 7 x 10^ km. A 



relationship for the radius of the companion is given by iKingI (jl988[ ) 



R, ~ 10^ km (1 + Xf'^ ( — (5) 



M, 







where X is the hydrogen fraction. For X = 0.7 we find Rc ~ 8x 10'^ km, which is substantially 
smaller than the implied size of the eclipsing region. It seems likely that the eclipses are 
being caused by an extended region of gas surrounding the companion. 



4.4. PSR J1701-3006F 



PSR J1701-3006F (hereafter M62F) is a 2.29 ms binary MSP in a circular orbit. Like 
M62E, it has a low minimum mass companion (~ 22 Jupiter masses) and occupies a re- 
gion in Mc,min--Pb phase space typical of non-eclipsing black widow pulsars. Indeed, despite 
full orbital coverage we see no evidence for flux variability as a function of orbital phase. 
Nonetheless, given the very low mass limit and highly circularized orbit, we still favor clas- 
sifying M62F as a black widow, an d the lack of eclipses probably indicate that the system is 
not being viewed close to edge-on (jFreirdl2005l ). 



4.5. PSR J1807-2459A 



PSR J1 807-2459A (NGC 65 44A) has already been discussed extensively by lD'Amico et al. 



teOOlal ) and lRansom et al.l ( 120011 ). It is a 3.06 ms MSP in a black widow system with a highly 
circular orbit and Mc^min = 0.009 Mq. We were able to reliably phase connect the data col- 
lected in 2009-2010 to data taken in 2011 as part of our Shapiro delay observations of PSR 
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J1807— 2500B (see below), as well as to older GBT data from 200#i. The orbital parameters 
of this system were already well determined by previous authors, however no phase- coherent 
timing solution had ever been derived or published. Our new timing solution includes, for 
the first time, a precise measurement of the position of the system and its spin-down. The 
orbit in our solution is in good agreement with previous analyses, but significantly more 
precise; it allowed us to measured Pb- Like M62B, the GR contribution to Pb is two orders 
of magnitude smaller than observed for a 1.4 Mq pulsar and minimum mass companion. 
However, there are no eclipses observed for NGC 6544A, so it is possible that the orbital 
inclination is significantly smaller than 90°, which would imply a higher companion mass 
and larger GR contribution. Nonetheless, for GR to contribute even 10% of the measured 
Pb, the companion would have to have a mass of ~ 0.06 Mq, corresponding to i = 9°. For 
a distribution of inclination angles that is fiat in cos i, there is only a ~ 1% probability of 
having i < 9°. It thus seems that we can safely rule out a significant contribution from 
GR to Pb. We can also rule out significant contamination from acceleration in the cluster 
because the observed Pb/Pb is several times larger than the maximum cluster acceleration. 
It thus seems lik ely that Pb can be explained by normal long-term black widow behavior 



( iNice et al.ll2000l ). NGC 6544A is offset from the center of the cluster by 4.6", or roughly 



1.5 core radii. At the distance of NGC 6544, this is only 0.06 pc. 

The observed P < 0, which, if it were intrinsic to the pulsar, would imply that the 
pulsar is spinning up. In reality, P is contaminated by the acceleration of the pulsar within 
the gravitational field of the cluster and Galaxy, and the fact that P < gives unambiguous 
evidence that the pulsar is on the far side of the cluster and accelerating towards Earth. 
This provides us with a probe o f the mass enc losed at the projected posi t ion of the pulsar 



and the mass-to-light ratio (T) (jPhinneylll993l ). Following iD'Amico et al.l ( 12002| ). 



Tv > 1.96 X 10^^ ^ ( ^^^^^ ) " (6) 

where Sv(< 0±) is the mean surface brightness interior the position of the pulsar. To 
calculate Sv(< 0±), we assume a constant surface brightness in the core of the cluster. 
The cluster acceleration, ac, is obtained by re-arranging Eq. [H In this case, the intrinsic 
spin-down of the pulsar is estimated by using the formula for characteristic age (tc), 

^ - — (7) 



P 2r, 



^*The difference between the pre- fit residuals of the 2004 data and the predictions of our 2009-2011 based 
timing solution are only 6% of pulse phase. Therefore, we are confident that we can phase connect all our 
data. 
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assuming a pulsar age of 10 Gyr. We find Pmt/P = 1-6 x 10^^^ s^^. We were unable to 
find a proper motion for NGC 6544 in the literature, but our results are not very sensitive 
to this — for example, if the cluster had a transverse velocity of 100 km s~^, it would affect 
our results at the 10% level. Combining all of this information, we find Ty > 0.072 Mq/Lq. 
Massive stellar remnants in the cores of GCs should give rise to a higher Ty (typically 
~ 3-4 M0/L0), so this result is very unconstraining. 



4.6. PSR J1807-2500B 



PSR J1807-2500B (NGC 6544B: IChandlei]l2003[ ) is the most intriguing system in our 
sample. As with NGC 6544A, were able to phase connect data spanning 2004-20110. The 
pulsar has a period of 4.19 ms and is in a highly eccentric orbit, with e = 0.747 (see Fig. 
[3]). The very high eccentricity of this system has allowed us to measure a very significant 
precession of periastron, u = 0.018319(12)° yr~^. We have also measured Shapiro delay in 
NGC 6544BP^. We discuss both measurements in more detail below. 



4.6.1. Possible Contributions to u 

The observed orbital precession could be due to any combination of tidal deformation of 
the companion, spin-orbit coupling, or GR effects, but for various reasons discussed below, 
we believe that it is due almost entirely to GR. 

Tidal deformation of the companion would require a MS or giant companion. This 
already seems unlikely based solely on the mass function — the minimum companion mass 
for Mp = 1.4 Mq is ~ 1.2 Mq, which is well above the turn-off mass in GCs (~ 0.8 Mq). 
Even if we use a smaller pulsar mass (say 1.2 Mq), the minimum companion mass is still 
~ 1.1 Mq. This rules out a MS companion, unless it was a blue straggler, but given how rare 
these are, this seems exceedingly unlikely. Any giants would have to be close to the turn-off 
mass and would experience Roche lobe overflow. The resulting gas in the system would give 
rise to eclipses, especially if the orbit is highly inclined (and as we show below, it is) and near 
conjunction. As it turns out, we see no evidence for eclipses at any orbital phase, including 
conjunction, so a giant seems unlikely. A blue straggler would also probably cause eclipses 



^^The difference between the pre- fit residuals of the 2004 data and the predictions of our 2009-2011 based 
timing solution are < 3% of pulse phase. Therefore, we are confident that we can phase connect all our data. 

^®To the best of our knowledge, this is the first detection of Shapiro delay in a GC pulsar. 
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or some other timing irregularities in the pulsar. In fact, NGC 6544B times remarkably well, 
with no unmodeled trends in the data and no need for an error factor to bring Xrcd — ^■ 
Hence, we can safely rule out a significant contribution to u by tidal deformation. 

Spin orbit coupling is a possibility if the companion is rapidly rotating, but it scales as 
\x/x\ times a geometrical factor, where x is the projected serni- maj or axis; the geometrical 



factor is expected to be < 10 in 80% of cases (e.g. iFreire et al.ll2008l . and references therein). 
Our current best limit on x implies that precession due to spin-orbit coupling should be 
smaller than our extremely small measurement uncertainties in uj 

For the reasons outlined above, we feel confident that the observed u is due almost 
entirely to GR. In this case the total mass of the system is given by 



5/2 



where Tq = GMqC ^, and Pb is the orbital period. Using our measured values, we find 
Mtot = 2.56763(59) Mq. 



4-6.2. Shapiro Delay Measurement 

The total system mass of NGC 6544B is fairly high for a pulsar binary system. When 
combined with the observed mass function, it implies a high orbital inclination and massive 
companion, making NGC 6544B an excellent candidate for a measurement of Shapiro delay. 
We observed the pulsar on ten occasions over a variety of orbital phases, including two eight 
hour tracks at or near orbital conjunctioij^. 



We used the DDH timing model developed by lFreire fc Wexl (120101 ) . which parameterizes 
the Shapiro delay as 

^ ^ YT7T=W 

where s = sini and r = TqMc are the traditional Shapiro delay parameters. In this param- 
eterization, quantifies the amplitude of the third harmonic and ^ the ratio of successive 



^^Because the orbital period is almost exactly ten days, conjunction moves slowly in sidereal time, and as 
such is not visible from a given telescope for long stretches. We were not able to schedule all our observations 
contiguously around conjunction due to this constraint. 
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harmonics of the Shapiro delay. These parameters are less correlated with each other and 
with other orbital parameters than are s and r, and provide a better description of the 
combinations of orbital inclination and companion mass allowed by the timing. From the 
measured values of and <^ we derive sinz = 0.99715(20) and = 1.02(17) Mq. On its 
own, the Shapiro delay does not provide precise measurements of the component masses, 
but the DDH model does not assume that u is relativistic, and thus does not make full use 
of the available information. 



To improve our mass estimates we used the DDGR model (ITaylor fc Weisberglll989l ). 
which assumes that GR correctly describes the system. In this model, Mtot and Mc are free 
parameters and all post-Keplerian (PK) parameters are derived from these measurements. 
From this we obtain M^f^^ = 2.57190(73) Mq, M°°^^ = 1.2064(20) Mq, and derive 
^DDGR _ 1.3555(21) Mq (see Table E] for the complete solution). This is the n aost p recise 
mass ever derived for an MSP, the previous being J1903+0327 (iFreire et all l201lh . To 
understand this exceptional precision and verify it, we created niaps from the DDH model, 
using the Keplerian parameters of the DDGR solution as a starting point. For each cosi and 
Mc, we calculated all five PK parameters {u, 7, Pb, ? and h^) assuming they are determined 
solely by GR; these were kept fixed while TEMP02 fit for all other parameters. The resulting 
values were recoded and the associated 2-D probability distributio n function (PDF) wa s 



calculated using the Bayesian analysis technique discussed in detail in ISplaver et al.l ( 120021 ). 
This 2-D PDF was then translated into the Mc-Mp plane using the mass function of our 
starting DDH solution. The original 2-D PDF was then collapsed onto the cos i and axes, 
with the derived 2-D PDF collapsed into the Mp axis (see red lines in Fig. [S]), generating 
1-D PDFs for the latter quantities. These are slightly asymmetric; they have medians 
and ±la percentiles at: cosi'^'"^ = Om7tofil, M^^^ = 1.3649l|];[^[J^^ Mq, and M^^'^^ = 
1.2068tomll Mq. The maximum probability occurs at cosi"^"'' = 0.095, M^'" = 1.3654 Mq, 
and M™*^^ = 1.2063 Mq. These are in excellent agreement with the best-fit results of the 
DDGR model, and confirm the value and the small uncertainty of the masses. 

The regions of the cos i-Mp plane (and derived Mp-Mc plane) allowed by each measured 
PK parameter are depicted in Fig. [51 This allows us to understand the reason for such a 
precise mass: the constraint provided by u (the total mass Mtot, according to GR) cuts the 
cosi-Mc (or Mp-Mc) regions allowed by the Shapiro delay very sharply. 

The measured companion mass strengthens the case made in §4.6.1l that the companion 
must be a massive white dwarf or second neutron star. The companion is the most massive 
known around a fully recycled pulsar. If this i s a double neutron star system, it is only the 



second known in GCs, the first being M15C (lAnderson et al.lll990l : iJacoby et al.ll2006l ). It 



is likely that the pulsar was recycled by a different companion, which was then ejected in 
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an exchange interaction. Deep Hubble Space Telescope images of the cluster may be able to 
detect a white dwarf companion, while a non-detection of the companion would strengthen 
the case for a double neutron star system. 

We have measured three PK parameters {u, and h^). Assuming GR is correct, all 
three agree on the same region of the Mc-Mp plane. As such, GR passes this test, albeit 
at relatively low precision. Our current best measurement of a fourth PK parameter, the 
gravitational redshift, is 7 = 0.026(14) s. This is consistent with the GR prediction of 
0.014 s. We plan to continue monitoring this system long-term, and in a few years expect to 
have a more precise measurement of 7 that will allow for a second test of GR. 



4.7. PSR J1823-3021C 



PSR J1823-3021C (NGC 6624C) has P = 0.406 s, unusually slow among globular 
clu ster pulsa r s but, surprisingly, the second such slow pulsar in NGC 6624. It was discovered 
by IChandlerl ( 120031 ): we have now measured, for the first time, P = 2.25 x 10~^^ ss~^. Using 
Eq. [21 we find Pc,max = 1-7 x 10~^^ s s~^, over an order of magnitude smaller than the 
measured P. The contributions from the Galaxy and the Shklovskii effect are even smaller 
still, so it is clear that the observed P is due almost entirely to the intrinsic pulsar spin- 
down. The implied characteristic age and surface magnetic field is ~ 2.8 x 10^ yr and 
3.1 X 10^1 G. This makes NGC 6624C, like NGC 6624B, similar to the "normal", non- 
recycled pulsars (NRPs) usually seen in the Galactic disk. As explained in ^ NRPs are 
typically assumed to form in core collapse supernova, which require massive stars that have 
not existed in GCs for billions of years. Since NRPs have lifetimes <C 10^ yr, core collapse 
supernovae cannot explain the presence of NGC 6624C and pulsars like it. The leading 
alternative explanation is electron capture supernovae (ECS). The kick velocities (fkick) that 
pulsars receive whe n they forrn via ECS are not well known, though there is evidence that 
vv^rk - 10 kms~^ flPfahl et al.l[2002[ iKitaura et aD bood : lOessart et aP l2006[ iMartin et aL 



2009 : Wong et al. 201o[ ). The escape velocity from the center of NGC 6624 is ~ 35 km s 



which effectively places an upper bound on the f kick that NGC 6624C received. This supports 
the notion that ECS kicks are much smaller than those induced by core collapse supernovae. 



GC NRPs can also ena ble useful statistical constraints on the properties of ECS (IBoyles et al 



2OIII : iLynch et al.llin prep.l ). NGC 6624C is one of only four GC NRPs known, making it an 
important addition to this rare family of pulsars. 
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4.8. Pulsars Without Full Timing Solutions 

As part of this work, wc discovered three previously unknown MSPs, J 1823-302 ID, E, 
and F (hereafter NGC 6624D, E, and F). However, we were unable to obtain full timing 
solutions for them. This was due to the very low signal-to-noise ratio of the detections of all 
three pulsars, due in part to their low flux densities, but also due to a data processing error on 
our part. As is typical, we de-dispersed the raw data and combined many frequency channels 
to reduce data volume (i.e. sub-banding), discarding the raw data afterward. However, we 
accidentally de-dispersed at the wrong DM for many of our observations, thereby adding 
significant dispersive smearing and degrading the signal from these already weak pulsars 
even further. Without the raw data, we are unable to remedy this error. PSR J1823— 3021A 
was bright enough that it could still be detected, and the dispersive smearing was not large 
enough to significantly impact the two long-period pulsars. Despite this, we could still 
determine some of the basic properties of NGC 6624D, E, and F. 

We are able to confirm that NGC 6624D and E are isolated MSPs with P = 3.02, and 
4.39 ms, respectively. NGC 6624D was detected on many occasions, but the TOA errors were 
quite large. We were able to measure the position with some accuracy, although without 
a reliable measurement of P, we are unable to break all the covariances in the fit, so these 
coordinates should be used cautiously. NGC 6624E was only detected in three observations 
but there was no variation in the apparent period of the pulsar between these observations, 
nor was there any evidence of acceleration. Both of these effects are expected if the pulsar 
were in a binary system. We were unable to phase connect any of our observations due 
to large relative errors on the measured period of the pulsar, and we cannot constrain its 
position or P. 

NGC 6624F is an eclipsing binary MSP with P = 4.85 ms. Like NGC 6624E, NGC 6624F 
was only detected three times. In each case, we were able to clearly measure variations in 
P, P and P arising from orbital motion. As it turns out, the apparent period of the pulsar 
was similar and the acceleration was positive in all three observations, indicating that the 
observations occurred at overlapping orbital phases near inferior conjunction. Because the 
pulsar had the exact same apparent period at some point during each of our observations, 
there must have been an integer number of orbits between these times. We calculated the 
exact moment during each observation when the pulsar had a reference barycentric period 
of precisely 4.85 ms through a Taylor expansion of P, P, and P and a bisection method. An 
orbital period of 0.8827 days implies exactly 106.0003 and 146.9996 complete orbits between 
our three observations. The fact that we were able to fit such a precise integer number of 
orbits gives us confidence that this is the true orbital period of the system. If NGC 6624F 
is in a circular orbit, each observation would mark points on an ellipse in the P-P plane 
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( iFreire et al.ll200ll ). With three observations and a good determination of the spin and or- 
bital period, we can, in principal, uniquely constrain the equation of this ellipse, which in 
turn can provide the intrinsic P, P]^, and a sin i/c. We attempted to perform this fit to our 
data using a least squares minimization routine and 0.01 It — s < asinz/c < 1000 It — s. We 
found a best-fit a sin i/c ~ 4.4 It — s, but this solution has poor predictive power (i.e., when 
we folded our data using this solution, the pulsar was not detected). We were also unable to 
construct a coherent timing solution in TEMP02 using this starting orbital solution. We be- 
lieve that the lack of coverage at multiple orbital phases is limiting our ability to accurately 
determine a sin i/c. It is also possible that NGC 6624F is not in a circular orbit. 



5. Conclusion 

We set out to obtain full timing solutions for nine previously known GC pulsars in 
M62, NGC 6544, and NGC 6624 (three of which have never been published). Our efforts 
were successful for six of these pulsars, and a seventh now has a partial solution. We 
have confirmed the nature of three black widow pulsars (M62E and F, and NGC 6544A). 
NGC 6544A has P < 0, providing unambiguous evidence that P is dominated by acceleration 
in the cluster potential, and providing us with a probe of Ty in the cluster core. We find 
Ty ^ 0.072 which, while not very constraining, is consistent with small amounts of low- 
luminosity matter such as massive black holes or other stellar remnants. We are also able to 
confirm that NGC 6624C belongs to a rare class of non-recycled GC pulsars similar to the 
slow pulsars found in the Galactic disk. 

The highlight of our work is unquestionably NGC 6544B. This binary MSP is in a 
highly eccentric orbit and exhibits clear orbital precession and Shapiro delay. Under the 
well justified assumption that this is due almost entirely to GR, we are able to obtain a 
most-probable total system mass of 2.57190(73) Mq, companion mass of 1.2064(20) M©, 
and pulsar mass of 1.3655(21) Mq. This is the highest mass companion known to orbit 
a fully recycled MSP, and raises the possibility that NGC 6544B is part of either double 
neutron star system, which could form via exchange interactions in the GC. 

We are deeply grateful to Andrea Possenti and Alessandro Corongiu for providing us 
with data collected by the Parkes Telescope that was used to check for phase connection in 
our solutions for M62D, E, and F. We would also like to thank an anonymous referee for 
helpful comments. R. Lynch acknowledges support from the GBT Student Support program 
and the National Science Foundation grant AST-0907967 during the course of this work. 
The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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Table 1. Summary of Pulsars in this Study 



Pulsar Previous Full New Full Reference 

Timing Solution? Timing Solution? 



M62 



J 170 1-3006 A 


Yes 




Yes 


A,B 


J1701-3006B 


Yes 




Yes 


B 


J1701-3006C 


Yes 




Yes 


B 


J1701-3006D 


No 




Yes 


C 


J1701-3006E 


No 




Yes 


C 


J1701-3006F 


No 




Yes 


C 


NGC 6544 


J1807-2459A 


No 




Yes 


A,D 


J1807-2500B 


No 




Yes 


C 


NGC 6624 


J1823-3021A 


Yes 




Yes 


E 


J1823-3021B 


Yes 




Yes 


E 


J1823-3021C 


No 




Yes 


C 


J1823-3021D 


No 




Yes 


Unpublished 


J1823-3021E 


No 




No 


Unpublished 


J1823-3021F 


No 




No 


Unpublished 



References. 



A-D'A 



C- IChandleij (l2003h : D- 



mico et all (l2001ah: B 



Ransom et al. 



Possenti et al. 


(20031 


— Biggs et al. 


1994) 



Table 2. Globular Cluster Properties 



ID 








„ d 
' c 




f 




(deg) 


(deg) 


(kpc) 


(arcmin) 


(mag arcsec"^) 


(kms-1) 


M62 


353.57 


7.32 


6.7 


0.22 


15.10 


13.3 


NGC 6544 


5.84 


-2.20 


3.0 


0.05 


16.31 


5.89 


NGC 6624 


2.78 


-7.91 


7.8 


0.06 


15.32 


6.0 



^Galactic longitude 
^Galactic latitude 
'^Distance 
'^Core radius 

"Central V-band surface brightness 
^Central velocity dispersion 



Note. — All properties are taken from 



av in M62 and NGC 6544 (IWebbinklll985l ). and NGC 6624 (jValenti et al.ll201l[ ). 



HarrisI |l996l 201 edition), except fo r 



Table 3. Parameters of the Previously Timed Pulsars in M62 



Parameter 


J1701-3006A 


J1701-3006B 


J1701-3006C 




Timing Parameters 


Right Ascension (J2000) 


17:01:12.50942(49) 


17:01:12.6677(10) 


17:01:12.86592(35) 


Declination (J2000) 


-30:06:30.173(36) 


-30:06:49.044(76) 


-30:06:59.415(21) 


Spin Period (s) 


0.0052415662037958(29) 


0.0035938521270515(34) 


0.0076128487111379(68) 


Period Derivative (ss~^) 


-1.301(15) X 10-1'' 


-3.483(20) X 10-19 


-6.413(10) X 10-20 


Dispersion Measure (pc cm~^) 


114.9654(84) 


115.21(38) 


114.5619(69) 


Reference Epoch (MJD) 


55038.0 


55038.0 


55038.0 


Span of Timing Data (MJD) 


53150-55200 


53315-55200 


53315-55200 


Number of TOAs 


80 


81 


74 


RMS Timing Residual (fis) 


6.4 


11.3 


15.1 


Error Factor 


1.6 


1.5 


1.2 


Binary Parameters 


Orbital Period (days) 


3.8059483732(58) 


0.14454541718(58) 


0.21500007119(51) 


Projected Semi-major Axis (It-s) 


3.4837397(19) 


0.2527565(28) 


0.1928816(33) 


Epoch of Periastron Passage (MJD) 


55040.03822060(36) 


55038.62576448(54) 


55037.92664751(74) 


Orbital Eccentricity 


< 1.8 X 10-8 


< 4.5 X 10-5 


< 7.9 X 10-8 


Longitude of Periastron Passage (deg) . . . 








Rate of Change of Orbital Period (lO-^^) 




-5.12(62)xl0-i2 




Derived Parameters 


Mass Function (M0) 


0.0031339643(51) 


0.000829814(27) 


0.0001666778(87) 


Minimum Companion Mass (Mq) 


0.2 


0.12 


0.071 


Offset from Cluster Center (core radii) . . . 


1.50 


0.138 


0.773 


Intrinsic Spin-down (ss^-"-) 


< 3.13 X 10-19 


< 4.4 X 10-20 


< 3.78 X 10-19 


Surface Magnetic Field (10^ Gauss) 


< 1.3 


< 0.4 


< 1-7 


Spin-down Luminosity (10^* ergs~^) 


< 8.6 


< 3.7 


< 3.4 


Chara<;teristic Age (10^ yr) 


> 0.27 


> 1.3 


> 0.32 


Mean 2 GHz Flux Density (/xJy) 


117 


78 


72 



Note. — Numbers in parentheses represent l-cr uncertainties in the last digits as determined by TEMP02, scaled such that 
the Xred ~ ^- timing solutions use the DE405 Solar System Ephemeris and the TT(BIPM2011) time standard. All values 
are reported in Bajrycentric Dynamical Time (TDB) units. Mean flux density estimates have a typical 10%— 20% relative error. 
Minimum companion masses were calculated assuming a 1.4 Mq pulsar. 



Table 4. Parameters of the Newly Timed Pulsars in M62 



Pcirameter 


J1701-3006D 


J1701-3006E 


J1701-3006F 




Timing Parameters 


Right Ascension (J2000) 


17:01:13.5631(15) 


17:01:13.2734(17) 


17:01:12.8222(15) 


Declination (J2000) 


-30:06:42.559(79) 


-30:06:46.89(11) 


-30:06:51.82(10) 


Spin Period (s) 


0.0034177704450548(17) 


0.0032337373331158(65) 


0.0022947270806719(29) 


Period Derivative (ss^^) 


1.257(24) X 10-19 


3.103(27) X 10-19 


2.221(18) X 10-19 


Dispersion Measure (pc cm~^) .... 


114.224(11) 


113.792(22) 


113.291(34) 


Reference Epoch (MJD) 


55038.0 


55038.0 


55038.0 


Span of Timing Data (MJD) 


54876-55200 


54876-55200 


54876-55200 


Number of TOAs 


45 


40 


47 


RMS Timing Residual (/us) 


6.5 


13.1 


19.5 


Error Factor 


1.1 


1.4 


1.3 


Binary Parameters 


Orbital Period (days) 


1.1179034034(82) 


0.1584774951(56) 


0.2054870422(72) 


Projected Semi-major Axis (It-s) . . 


0.9880185(50) 


0.0701589(55) 


0.0573400(52) 


Epoch of Periastron Passage (MJD) 


55037.88706413(93)=- 


55037.9364630(50) 


55038.0892403(49) 


Orbital Eccentricity 


0.0004122(67) 


< 1.9 X 10-4 


< 3.4 X 10-" 


LoiigitiuU* (if I\'riasl roil l^issa.y;t* (dv 


k) 101. fid, 3) 






Derived Parameters 


Mass Function (Mq) 


0.000828647(13) 


1.47637(35) X 10-5 


4.7939(13) X 10-8 


Minimum Companion Mass (Mq) . 


0.12 


0.031 


0.021 


Offset from Cluster Center (core radii) 0.912 


0.503 


0.185 


Intrinsic Spin-down (s s~^) 


< 3.4 X lO-i'' 


< 5.56 X 10-19 


< 4.14 X 10-19 


Surface Magnetic Field (10^ Gauss) 


< 1.1 


< 1-4 


< 0.99 


Spin-down Luminosity (10''^ ergs^-'- 


) ■ < 34 


< 65 


< 140 


Characteristic Age (10'* yr) 


> 0.16 


> 0.092 


> 0.088 


Mean 2 GHz Flux Density (/iJy) . . . 


69 


39 


52 



^For M62D we give the epoch of the ascending node (Tasc), which is more precisely measured in the ELLl timing model. 
Note. — See the notes to TabIe|3]for more details on these timing solutions. 



to 
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Table 5. Parameters of NGC 6544 A 



Parameter J1807-2459A 



Timing Parameters 

Right Ascension (J2000) 18:07:20.355604(16) 

Declination (J2000) -24:59:52.9015(65) 

Spin Period (s) 0.003059448798020229(30) 

Period Derivative (ss"!) -4.3352(26) X lO'^l 

Dispersion Measure (pc cm"^) 134.00401(58) 

Reference Epoch (MJD) 55243.0 

Span of Timing Data (MJD) 53315-55760 

Number of TOAs 296 

RMS Timing Residual (/xs) 1.4 

Error Factor 1.4 



Binary Parameters 

Orbital Period (days) 0.071091483516(29) 

Projected Semi-major Axis (It-s) 0.01222393(12) 

Epoch of Periastron Passage (MJD) 55242.99436660(24) 

Orbital Eccentricity < 1.1 X 10"'' 

Longitude of Periastron Passage (deg) ... 

Rate of Change of Orbital Period (lO'^^) -1.142(62) X IQ-^^ 



Derived Parameters 

Mass Function (M©) 3.88043(11) X 10"^ 

Minimum Companion Mass (Mq) 0.0092 

Offset from Cluster Center (core radii) ... 1.31 

Intrinsic Spin-down (ss-i) < 2.58 X 10~i^ 

Surface Magnetic Field (10^ Gauss) < 0.9 

Spin-down Luminosity (10'^'' ergs~^) .... < 36 

Characteristic Age (10^ yr) > 0.19 

Mean 2 GHz Flux Density (/.tjy) 605 

Rotation Measure (radm~-^) 160.4(6) 



Note. — See the notes to Table [3l for more details on this timing solution. 
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Table 6. Parameters of NGC 6544B using DDGR Model 



Parameter 


J1807-2500B 




Timing Parameters 


Right Ascension (J2000) 


18:07:20.871209(53) 


Declination (J2000) 


-25:00:1.915(17) 


Spin Period (s) 


0.00418617720284089(25) 


Period Derivative (ss~^) 


8.23245(18) x 10"^" 


Dispersion Measure (pc cm~'^) 


137.1535(20) 


Reference Epoch (MJD) 


54881.352079 


Span of Timing Data (MJD) 


53315-55760 


Number of TOAs 


278 


RMS Timing Residual (/xs) 


5.4 


Error Factor 


1.0 


Binary Parameters 


Orbital Period (days) 


9.9566681588(27) 


Projected Semi-major Axis (It-s) 


28.920391(44) 


Epoch of Periastron Passage (MJD) 


54881.34735776(22) 


Orbital Eccentricity 


0.747033198(40) 


Longitude of Periastron Passage (deg) . . . 


11.334600(18) 


Total System Mass (Mq) 


2.57190(73) 


Companion Mass (Mq) 


1.2064(20) 


Derived Parameters 


Pulsar Mass (Mq) 


1.3655(21) 


sin i 


0.9956(14) 


Mass Function (Mq) 


0.2619794(12) 


Rate of Change of Orbital Period (10" 


-0.027921(62) 


Rate of Periastron Advance (degyr~^) . . . 


0.0183389(35) 


Gravitational Redshift (s) 


0.014418(40) 


Offset from Cluster Center (core radii) . . . 


4.05 


Intrinsic Spin-down (ss~^) 


< 2.25 X 10-1^ 


Surface Magnetic Field (10^ Gauss) 


< 0.98 


Spin-down Luminosity (10'^* erg s~^) .... 


< 12.0 


Characteristic Age (10^ yr) 


> 0.29 


Mean 2 GHz Flux Density (/xjy) 


88 


Rotation Measure (rad m~'^) 


157(1) 



Note. — In the DDGR model, Mtot and Mc are the only free PK pa- 
rameters. See the notes to Table |3] for more details on this timing solution. 



Table 7. Parameters of the Pulsars in NGC 6624 



Parameter 



J1823-3021A 



J1823-3021B 



J1823-3021C 



J1823-3021D 



Timing Parameters 



Right Ascension (J2000) 18:23:40.48701(36) 

Declination (J2000) -30:21:40.127(44) 

Spin Period (s) 0.0054400041632727(31) 

Period Derivative {ss'^) 3.3750(20) X lO'^** 

Dispersion Measure (pc cm""^) 86.8797(56) 

Reference Epoch (MJD) 55049.0 

Span of Timing Data (MJD) 54898-55200 

Number of TOAs 33 

RMS Timing Residual (/xs) 4.6 

Error Factor 1.1 

Derived Parameters 

Offset from Cluster Center (core radii) 0.115 

Intrinsic Spin-down (s s'^) < 3.64 X 10~^^ 

Surface Magnetic Field (10^ Gauss) . . < 45 

Spin-down Luminosity (10'^'* ergs~^) . < 8900 

Characteristic Age (10^ yr) > 0.00024 

Mean 2 GHz Flux Density (^Jy) 79 



18:23:41.5455(23) 
-30:21:40.94(47) 
0.3785964921974(27) 
2.99(14) X 10-^^ 
87.01(10) 
55049.0 
54898-55200 
33 
55.8 
1.0 



18:23:41.1516(36) 
-30:21:38.45(80) 
0.4059359629899(53) 
2.240(25) X 10-1'' 
86.88(24) 
55049.0 
54898-55200 
31 
104.9 
1.1 



18:23:40.5312(71) 
-30:21:43.66(35) 
0.003020060260978(59) 

86.8(1.1) 
55050.0 
54898-55200 
10 
90.0 
3.9 



3.74 
< 3.45 X 10- 
< 120 
< 0.0025 
> 0.17 
76 



2.24 
< 2.33 X 10-1® 
< 310 
< 0.014 
> 0.028 
37 



2.33 



29 



Note. — 



See the notes to Table |3] for more details on these timing solutions. 
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Fig. 1. — Calibrated pulse profiles for NGC 6544A and B showing the integrated pulse 
profile (top solid lines) and the fraction of linear (second solid line) and circular (dotted line) 
polarization. The top panels show the polarization angle. The profiles have been rotated by 
an arbitrary amount. These data are from our Shapiro delay observations for NGC 6544B 
(see §4.6p and were taken at 1.4 GHz in a coherent de-dispersion mode. The sampling time 
was 10.24 fis (0.3% and 0.2% of pulse phase for NGC 6544A and B, respectively). 
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Fig. 2. — Average 2 GHz pulse profiles for all pulsars in these clusters. Binary pulsars are 
denoted by an asterisk. These profiles formed the basis of the standards used in our timing 
analysis, and were created by summing all the detections of a given pulsar. In the case of the 
pulsars in NGC 6624, only data de-dispersed at the correct DM was used for making these 
plots (see text for more details). The bars indicate the contribution of dispersive smearing 
to the profile width. 



-31 - 



0.349 



-0.116 



0.497 



0.104 



0.035 



-0.035 



0.038 



-0.013 



M62A (Pnrb = 3.8 days) 




M62C (Porb = 5.2 hr) 




M62E (Porb = 3.8 hr) 




NQC 6544A (P^b = 1.7 hr) 
1 



0.457 



0.152 



-0.152 



-0.457 



0.220 



0.073 



- -0.073 



-0.220 



0.047 



0.016 



-0.016 



-0.047 



2,304 



M62B (Pcrb = 3.5 hr) 




M62D (Porb = 1.1 days) 




M62F (Porb = 4.9 hr) 




NGC 6544B (Po^b = 10.0 days) 




0.5 1.0 1.6 

Mean Anamoly 



2.0 



Fig. 3. — The change in observed rotational period due to the Doppler effect as a function of 
mean anomaly. The circles show the observed periods. For clarity, the data are duphcated 
from orbital phase 1-2. Note that the high eccentricity of NGC 6544B causes the large 
deviation from a sinusoidal modulation of the period. 
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Fig. 4. — Post-fit timing residuals for the pulsars in these clusters. Only unambiguously 
phase connected TOAs are displayed. The TOAs around MJD 55600 for NGC 6544A and 
B are from recent Shapiro delay observations of NGC 6544B. 
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Fig. 5. — Constraints on the inclination angle and masses in NGC 6544B from maps and 
the DDH timing model. Our measurement of uj implies a total mass of 2.57190(73) Mq, 
indicated by the black lines. The 1-a uncertainty is so small that it is invisible at this 
scale. The solid purple lines indicate the constraints from and the dashed purple lines 
indicate the constraints from q. Left: cosi-Mc plot. The gray region is excluded by the 
condition Mp > 0. Right Mp-Mc plot. The gray region is excluded by the condition 
I sin^l < 1. The black contours include 68.3% of the total probability of a 2-D probability 
distribution function (PDF), calculated from a map that used only the Shapiro delay 
to constrain the masses. The red contour levels further assume that GR is correct and 
is the sole cause of u. Top and right marginal plots: 1-D PDFs for cosi, Mp and M^, 
obtained by marginalization of both 2-D PDFs. The second PDF (in red) yields median 



masses of M™^'^ 



at M^'^^ 



1.3654 M0 and 



+0.0017 y^^^ ^mcd 

1.2063 M0. 



1.3649_o,oo22 ^^^05 



1.2068 



+0.0022 
-0.0016 



Mq. The probability peaks 
These are in excellent agreement with the 



masses obtained from the DDGR timing solution, given in Table El 



